Magnetic polarons in doped ID antiferromagnetic chain 
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The structure of magnetic polarons (ferrons) is studied for an ID antiferromagnetic chain doped 
by non-magnetic donor impurities. The conduction electrons are assumed to be bound by the 
impurities. Such a chain can be described as a set of ferrons at the antiferromagnetic background. 
We found that two types of ferrons can exist in the system. The ground state of the chain corresponds 
to the ferrons with the sizes of the order of the localization length of the electron near the impurity. 
The ferrons of the second type produce a more extended distortion of spins in the chain. They are 
stable within a finite domain of the system parameters and can be treated as excitations above the 
ground state. The ferrons in the excited states can appear in pairs only. The energy of the excited 
states decreases with the growth in density of impurities. This can be interpreted as a manifestation 
of an attractive interaction between ferrons. 
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I. INTRODUCTION 

It is now commonly accepted that the tendency to- 
ward phase separation is of fundamental importance for 
the physics of manganites and related compounds^^. 
The self-trapping of charge carriers is the most widely 
discussed type of phase separation, first predicted in 
the seminal paper of Nagaevi In such phase-separated 
states, charge carriers are confined within small ferro- 
magnetic metallic droplets (magnetic polarons or ferrons) 
located in an insulating antiferromagnetic matrix. It is 
usually assumed that the region of the perturbed spins 
around a ferron is rather narrow that is confirmed by nu- 
merical calculations for some models^. However, as first 
pointed out by De Gennes^, the distortion of the mag- 
netic order around a magnetic defect (e.g. a ferron) may 
decay slowly with the distance. The possible existence 
of ferrons with such extended spin distortions was an- 
alyzed for one-dimensional (ID) antiferromagnetic chain 
in Refs.0EJ These calculations performed for an isolated 
impurity show that the characteristic length of the dis- 
torted spin surrounding can be much larger than the size 
of the trapping region. 

The problem arises what kind of ferrons should be ex- 
pected at different values of the parameters of the system. 
To study this problem, we consider in this paper ID anti- 
ferromagnetic chain doped with non-magnetic donor im- 
purities of finite density. We find that there are two pos- 
sible magnetic structures: one characterized by extended 
spin distortions and another one with a narrow distorted 
region. For the relationship between parameters char- 
acteristic of manganites, the ground state for a regular 
distribution of impurities corresponds to the ferrons with- 
out extended spin distortions. The ferrons with extended 
spin distortions can appear as excited metastable states. 
We found that due to the overlapping of extended spin 



distortions, an attractive interaction among these ferrons 
arises in the system. This interaction could favor a phase- 
separation process, that would lead the system to an in- 
homogeneous state with clusters of ferrons, as observed 
in the experiment. 



II. THE MODEL 

In this section, we find the magnetic structure of the 
ground state of an ID antiferromagnetic chain doped 
by non-magnetic donor impurities. We consider an 
one-dimensional chain of antiferromagnetically coupled 
local spins, the y-axis being the axis of the chain. 
We do not considered here specific effects related to 
one-dimensionality and start from the antiferromagnetic 
structure characteristic of the mean-field approximation. 
Non-magnetic impurities occupy the sites of the chain 
at half-integer positions. The conduction electrons are 
assumed to be in a bound state of the impurity electro- 
static potential. It is also assumed, for simplicity, that 
the wave function of the conduction electron extends only 
over the two neighbouring magnetic ions near an impu- 
rity. Further on, we assume that for a given density of 
conduction electrons n, the impurities are periodically 
distributed along the chain, and L = 1/n > 2 being the 
distance between neighboring impurities (in units of the 
lattice constant). We consider the following Hamiltonian 
of the system 



H = H el + J' S g S g+1 - K' (S x g ) 2 , (1) 
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where 

H c \ = ~t^2 (a+ L a a iL+1 , a + aj L+l a a lL .^ - 

i,tj 

+ a iL+i, CT (f^iL+ij ^ ff) a lL+1:a , | (2) 

In Eqs. and (J2J, 5 9 is the spin of the magnetic ion lo- 
cated at site g, treated as classical vector, symbols a+ ff , 
a g _cr denote the creation and annihilation operators for 
the conduction electron with spin a at site g, and a are 
Pauli matrices. The second and third terms in Eq. (JTJ 
are the antiferromagnetic exchange between local spins 
and the magnetic anisotropy energy, respectively. The 
two terms in H e \ describe the kinetic energy of the con- 
duction electrons bounded by the impurities, which are 
located between sites with g — iL and g — iL + 1 (i 
is an integer), and the Hund-rule coupling between the 
conduction electrons and the localized spins. Parameters 
A, t, J', and K 1 of Hamiltonian (JTJ are considered to 
be positive and the x axis is the easy axis. The energy 
of Coulomb interaction between the conduction electron 
and the impurity is an additive constant and, for this 
reason, omitted in the calculation. 

Hamiltonian JIJ is applicable for the description both 
of wide-band and double-exchange (like manganites) 
magnetic semiconductors. The hierarchy of parameters 
in wide-band case is t ^> AS > J' while for the double- 
exchange case we have AS 3> t > J'. In both cases K' is 
assumed to be the smallest parameter, which corresponds 
in the most cases the experimental situation. 

It is well known that due to the Hund-rule coupling 
A Hamiltonian J2J) favors the deviations from the ideal 
antiferromagnetic arrangement of spins. The deviations 
are the most pronounced for the sites neighboring to im- 
purities. This can be illustrated by the diagonalization 
of Hamiltonian H c \ l|2"jl. which can be easily performed. 
As a result, we have that the lowest eigenvalue of H c \ 
corresponding to the ground state and the low-lying ex- 
citations can be written as 



E cl = -- Ja 2 S 2 + 4* 2 + iASt cos (|) , (3) 

i ~ 

where Vi is the angle between vectors So, and Sil+i- We 
can see that the lowest value of energy E e \ corresponds to 
the parallel orientation of spins nearest to the impurity. 
So, we get a "seed" ferron centered at the impurity. 

To the first approximation, both wide-band and 
double-exchange cases can be treated in the same man- 
ner. At t <C AS, the first-order term in E e \ with respect 
to t/AS has the form 

Be! , (4) 
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FIG. 1: The directions of local spins at sites with g = 0, 1, and 
g = 2, corresponding to the first (left panel) and second (right 
panel) configuration of local spins. The symmetry condition 
for the first (second) case is 9\ = 8q (8i = n — 0o). 



where N\ mp is the number of impurities. In the case 
t S> AS, the similar expression can be obtained by inter- 
changing t < — > AS/ 2. Further on, we discuss the double 
exchange case (@J, however, the results for the wide-band 
case can be obtained by the aforementioned change of the 
parameters. 

Let us consider now the configuration of the local spins 
corresponding to the ground state. Due to the rotational 
symmetry of the Hamiltonian with respect to the easy 
axis, all local spins in the ground state should lie in a 
plane containing the easy axis. The deviations of any 
local spin from this plane results in the growth of the en- 
ergy of the system. We assume without loss of generality 
that it is the xz plane. 

The direction of a local spin S g in the xz plane can 
be characterized by angle 9 g with respect to z axis 
(see Fig. [TJ. Then, the spins located between any 
two neighboring impurities can be represented as S g — 
S ((-l) 9 ~ s sin 9 g ,Q, cos 6 g \ where < g < n, and s = 

or 1. The passing through a ferron changes such a se- 
quence of spins to another one. Since the impurities are 
periodically arranged along the chain, the magnetic or- 
dering of local spins in the ground state has to be also 
periodic. There exist several ways to join the solutions 
at impurities, which preserve the periodicity. The corre- 
sponding relationships can be written as Sf£ = ±Sf£ +v 
However, only two of them give the energy gain and cor- 
respond to the ferron-type solutions. The first one is 
Sf L = —Sf L+1 , Sl L — S? L+1 . This means that the mag- 
netic moments of ferrons are directed along the z-axis. 
In the second case, the magnetic moments of ferrons are 
directed along the easy axis, and spins close to the ith im- 
purity satisfy the equality Sf L = S? L+1 , S z lL = -S? L+1 . 
The configuration of local spins corresponding to the first 

case is S g — S ((— I) 9 1 sin 6 g , 0, cos 9 g ^j for all g, and the 

symmetry condition is Oil — 9iL+i- The configuration of 
local spins corresponding to the second case can be writ- 
ten in the form S g — S ((— l) 9 1 1 sin6* g , 0, cos# ff ^ for 

iL + 1 < g < (i + 1)L, and the symmetry condition in 
this case is 9iL — tt ~ 9iL+i- Both configurations are 
schematically shown in Fig. ^ 

In the ID antiferromagnetic chain under study the 
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magnetic moments of neighboring ferrons are either par- 
allel or antiparallel for odd and even L, respectively. If L 
is odd integer, the periodicity condition is Q g +L = 0g, 
and the magnetic moments of all ferrons are parallel 
each other. For even L, the periodicity condition is 
9 g + l = 7r — Og and the magnetic moments of neighboring 
ferrons are antiparallel. 

Taking into account the periodicity of the problem 
and the aforementioned symmetry conditions, we readily 
come to the relationships between the canting angles 
in Eqs. and (0J and the angles On,, 9iL+i- For any L, 
these relationships can be written in the following form: 



cos 



fVi 

V 2 



cos 
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± 



'iL + l 



(5) 



where plus (minus) sign corresponds to the first (second) 
configuration. 

Now let us find the form of magnetic ordering for local 



spins and the energy of the system corresponding to each 
configuration. Substituting the expressions for S g and 
Eq. (J5J for canting angles i/i, into the Eqs. Q and 
we obtain the energy of the system in the following form 



E = J cos (6 g + 6 g+ i 

gjtiL g 

9iL ± OiL+1 



J2[Jcos(e lL ±e lL+1 )-2F 



cos 



,(6) 



where we use the notation J = J'S 2 , K = K'S 2 , and 
F = t/2 in the case t < AS, or F = AS /A for t > AS. 
The plus (minus) sign in Eq. (|SJ| corresponds to the first 
(second) configuration of local spins. The constant term 
in J3J is omitted here. 

To find the spin structure of the ground state, one has 
to minimize the energy © with respect to the angles 6 g . 
A set of nonlinear equations is obtained 



J sin (0 g + 6 g+1 ) + J sin (0 s _i + g ) + K sin (20 g ) = g ^ iL + 1, (i + 1) L 

J sin (8 lL+1 + 9 lL+2 ) ± J sin (8 lL ± 6 lL+1 ) + if sin (26 lL+1 ) = 2F-^-^ cos (^) 

J sin (%+i )L ± %+i)l+i) + J sm + %+i)l) + K sin (29 {l+1)L ) = 2F .^ d cos (^P 1 ) 

i 



Using the periodicity conditions 6*;^ = 0u + ul, 
0(i+i)L+i = ^x+i for odd or 6» iL = 7T - 0(i +1 )L, 
Q(i+i)L+i = Tt — 8iL+i for even L, the number of equations 
in system (J7J becomes equal to L. After this, the system 
of equations is solved numerically at different values of 
F/J and K/ J assuming that K/ J <C F j J. 1. As a result, 
we obtain the magnetic structure of the system of local 
spins with donor impurities. This structure is discussed 
in the next section. 



III. RESULTS 

The numerical analysis of the solutions of Eqs. (J7J 
demonstrates that two types of ferrons schematically 
shown in Fig. ^ can exist in the system under study. 
For relatively large F/J > 1.5 the ground state corre- 
sponds to the ferrons of the second type (see right panel 
in Fig.^). For F/J > 2, we have a classical ferron con- 
figuration 4 , namely, ideal ferromagnetic ordering inside 
the localization region of the doped electron (two lattice 
constants in our case) and ideal antiferromagnetism out- 
side this region. The solution corresponding to the first 
configuration of the local spins (see left panel in Fig. 
is characterized by a long-range deviation from the ideal 
antiferromagnetic ordering and can be considered as an 
excitation above the ground state. This spin structure is 



illustrated in Fig. [21 where we plot the directions of local 
spins between two impurities for n = 1/21. 

For smaller F/J, the ground state of the system corre- 
sponds to the first type of ferrons. However, in this limit 
the difference between the two types of ferrons is not so 
drastic. Both types of ferrons create a distortion of an- 
tiferromagnetic ordering outside the electron localization 
region. Moreover, the characteristic size of the distorted 
region is larger for the ferrons of the second type. There- 
fore, we come to the conclusion that for any value of the 
ratio F / J the ground state of the system corresponds to 
the ferrons with smaller radius of spin distortions. 

Let us consider in more detail the case F/J > 2 
(typical of manganites). Let us refer the ferrons with 
and without extended spin distortions as "coated" and 
"bare" , respectively. The system of "bare" ferrons forms 
the ground state, whereas, the system of "coated" ferrons 
can form a metastable state. In Figs. [3] and 01 we show 
the energy of the "coated" ferron and the angle of the 
local spin in the trapping region, 9q as functions of the 
inverse impurity density, l/n. The ferron energy i? po i is 
defined as the difference between energies of the system 
with and without doped electrons divided by the num- 
ber of impurities, that is, E vo \ — E/N imp + L(J + K). 
The energy of the "bare" ferron does not depend on the 
density of impurities (up to l/n = 1/2) and equals to 
Epal = —2(F — J). 
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FIG. 2: Magnetic ordering for the first configuration of 
local spins between two non-magnetic impurities (dots) for 
n = 1/2 f. This structure is repeated between any two non- 
magnetic impurities of the chain. It can be seen that in the 
region occupied by the conduction electrons (the first and 
the last magnetic sites), we have almost ferromagnetic order- 
ing, while in the remaining part of the chain, there exists a 
distorted antiferromagnetic ordering. The parameters of the 
Hamiltonian are F = 3, K = 2.5 1(T 2 (both in J units). 
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FIG. 3: Energy of "coated" ferron versus the inverse density 
of impurities L — 1/n. The values of parameters are F = 3, 
and K = 2.5 10 -2 (in J units). The energy of a "bare" ferron 
equals to E pol = -2(F — J) = -4. At L < L CI = l/n cr = 
12 (vertical line), the "coated" ferrons become unstable (see 
figures and text below). 



In Fig. |SJ we plot the energy of the "coated" ferron as 
a function of the canting angle v§ at different impurity 
densities. At small densities the minimum of the ferron 
energy corresponds to nonzero value of isq, whereas at 
the densities exceeding some critical value n cr the ferron 
energy has a minimum at isq = 0. 

The numerical analysis of Eqs. Q shows that for im- 
purity densities higher or equal than some critical density 
n CT , the angle #0 (which is one-half of the canting angle 
vq) for the "coated" ferrons is equal to zero (see Fig. 0}. 



0,25-, 



0,20- .* 

0,15- . 

0,10- 

0,05- 

• 

0,00- ••••••• 

4 8 12 16 20 24 28 32 36 40 

L=1/n 

FIG. 4: Angle 6q of the local spin and the z axis in the 
trapping region versus the inverse density of impurities, 1/n. 
This angle is one-half of the canting angle between spins in 
the trapping region. When the density of impurities is greater 
than or equal to the critical value n cr = 1/12, the canting 
angle is close to zero, and the state corresponding to "coated" 
ferrons becomes unstable. This is a result of the overlapping of 
the spin distortions caused by neighboring ferrons. The values 
of the parameters in the Hamiltonian are F = 3, K = 2.5 10 -2 
(in J units). 
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FIG. 5: Energy of the "coated" ferron versus canting an- 
gle at different impurity densities. When n > n CI the E po \ 
is a monotonic function of vo, and it has a minimum at the 
angle 1/0 = corresponding also to the system of "bare" fer- 
rons, which has a lower energy. This means that the system 
of "coated" ferrons is unstable, because there is no energy 
barrier between these two states. 

Moreover, all spins in the chain become perpendicular to 
the easy axis. Thus, at n > n cr , the periodic structure of 
"coated" ferrons becomes unstable. Indeed, such ferrons 
become structurally identical to "bare" ferrons but have 
higher energy. Any rotation of this spin structure toward 
the easy axis leads to the monotonic decrease of its en- 
ergy. In other words, at small densities each "coated" 
ferron creates a long-range spin distortion of the antifer- 
romagnetic background around it. This long-range spin 



5 



distortion stabilizes the "coated" ferrons. The radius of 
extended spin distortion created by a single "coated" fer- 
ron can be estimated as tq r* l/2n cr . When the density 
increases, the distortions created by neighboring ferrons 
start to overlap, the corresponding energy barrier lowers, 
and disappears at n = n CI . 

The critical density increases with K since any devi- 
ation of the spins from the easy axis become less favor- 
able with the growth of anisotropy energy. The value of 
n = n CI decreases with the increase of ratio Fj J since 
the smaller canting angles become more favorable with 
the growth of the electron kinetic energy t compared with 
the exchange integral J'. Note that for small values of 
F/J<1 both types of ferrons remain stable up to the 
density n = 1/2. 

Let us consider now the low density limit n — > 0. 
As can be seen from Figs. [3] and the ferron en- 
ergy and angle 9q are practically density independent 
at n <C n CI . Therefore, in zero approximation, we can 
consider "coated" ferrons as isolated objects. The nu- 
merical analysis shows that the difference between ener- 
gies of "coated" and "bare" ferrons AE, weakly depends 
on the ratio F/ J. Following Ref. 0, it can be shown in 
continuous approximation that AE ~ \/8KJ, which is 
confirmed by computational results. Since the ground 
state of the system corresponds to "bare" ferrons, the 
ferrons with extended spin distortions can be considered 
as elementary excitations with the energy AE. Note, 
that these excitations can appear only in pairs since a 
creation of a single "coated" ferron strongly disturbs the 
ground state and costs the energy at least of the order 
of J. It is easy to see that to create a pair of neighbor- 
ing "coated" ferrons, one needs to overcome the energy 
barrier of the order of K/n (because it is necessary to 
change the orientation of all spins between two ferrons 
with respect to easy axis). So, at low densities, this pair 
can be considered as a metastable state. 

The energy of a ferron in excited state AE is propor- 
tional to the square root of K and, for the case K <C J, 
it can be rather small in comparison with J. There- 
fore, even at low temperatures T -C T/v ~ J, one can 
expect to have a relatively large number of ferrons in 
excited state. Using conventional formulas for thermal 
averages and having in mind that "coated" ferrons arise 
in pairs, we have for the averaged density n\ of excited 
ferrons (ui/t = n/(l + cxp{AE/kT}). This approach is 
valid only when the density of impurities is less than ro cr , 
otherwise the ferrons could not be considered as weakly 
interacting particles. However, within the order of mag- 
nitude the estimates for (n±)T and AE are valid up to 
n ~ 2n cr . Then, at n > 2n cr the excitations under discus- 



sion become unstable if the temperature exceeds a certain 
critical value kT a ~ AE/\n(n/n cr - 1). At T > T cr the 
density of "coated" ferrons exceeds n cr . 

The "coated" ferrons can be treated as isolated objects 
only in the limit of n — ► 0. At finite n, it is necessary 
to take into account the density dependence of the fer- 
ron energy. Since the energy of the "coated" ferrons de- 
creases with the growth of their density, the interaction 
between them is attractive. This interaction arises due 
to the overlapping of the spin distortions of neighbor- 
ing ferrons and its characteristic length far exceeds the 
lattice constant. Although the calculation presented in 
Fig. |3 is valid only for a periodic distribution of the im- 
purities, it is natural to expect that the similar attractive 
interaction should take place for random distribution of 
"coated" ferrons. The detailed calculation of an attrac- 
tive interaction between ferrons in a general case is the 
subject of further research. We just note here that this 
attractive interaction can be, in principle, the possible 
mechanism for formation of a long-range phase-separated 
state involving large clusters of ferrons. 



IV. CONCLUSION 

We studied the magnetic structure for ID antiferro- 
magnetic chain doped by homogeneously distributed non- 
magnetic donor impurities. The ground state of the sys- 
tem corresponds to a set of ferrons with characteristic 
size of the order of the localization length of electrons 
near the impurity. The elementary excitations in the 
system are ferrons with extended long-range spin distor- 
tions. These excitations can be unstable for some range 
of parameters of the system, in particular, at sufficiently 
high impurity density. An attractive interaction exists 
between ferrons in excited state due to the overlapping 
of the extended spin distortions of neighboring ferrons. 
This attractive interaction can lead to a phase-separated 
state with large clusters of ferrons. 
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